Double-beta decay of 130 Te to the first 0+ excited state of 130 Xe with CUORICINO 
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The CUORICINO experiment was an array of 62 Te02 single-crystal bolometers with a total 130 Te 
mass of 11.3 kg. The experiment finished in 2008 after more than 3 years of active operating time. 
Searches for both 0v and 2u double-beta decay to the first excited + state in 130 Xe were performed 
by studying different coincidence scenarios. The analysis was based on data representing a total 
exposure of N( 130 Te)-t=9.5 x 10 25 y. No evidence for a signal was found. The resulting lower limits 
on the half lives are Tf ( 130 Te -> 130 Xe*) > 1.3 x 10 23 y (90% C.L.), and if ( 130 Te -^ 13 ° Xe*) > 

2 2 

9.4 x 10 23 y (90% C.L.). 

PACS numbers: 23.40.Hc, 23.40. Bw, 21.10.Bw, 27.60.+j 



I. INTRODUCTION 
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Two-neutrino double-beta (2vf3f3) decay and neutrino- 
less double-beta {Qvf3fi) decay have been known for over 
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70 years now [H, (a recent review can be found in 0]). 
While experimental evidence for 2i//3/3-decay has been 
found there is still no observation for the O^/3/3-decay, 
however several limits for the half-life have been set in 
the past with values greater than 10 21 y. In both of these 
processes the lifetime is proportional to the square of the 
Nuclear Matrix Elements (NME). Two neutrino double 
beta decay has been detected in ten nuclei on the ground 
state of the daughter nucleus and in two nuclei on the 
excited state of it, and the corresponding extracted val- 
ues for the NME are in reasonable agreement with the 
theoretical expectation. In the case of O^/3/3-decay their 
value is very important since it plays the same role in 
the prediction of the decay time as mpp, the effective 
neutrino mass 

The CUORICINO experiment was an array of 62 Te0 2 
bolometers operated at a temperature of about 10 mK. 
A bolometer 0, [1] detects an energy release as a tem- 
perature rise in the absorber crystal. Thermal pulses 
are converted into electric signals by means of neutron 
transmutation doped (NTD) thermistors @, which are 
coupled to each absorber. CUORICINO was organized 
in 13 planes. All of these planes were composed of four 
crystals with dimensions of 5 x 5 x 5 cm 3 and a mass of 
790 g each, except for the 11 th and 12 th (from top to bot- 
tom) . Each of these two particular planes had 9 crystals 
with dimensions of 3 x 3 x 6 cm 3 and a mass of 330 g. Two 
of these smaller crystals were enriched to 82.3% of 128 Te 
and two others to 75% of 130 Te. All the other crystals 
had the natural isotopic abundance of 130 Te (33.8%). A 
monthly calibration was performed using a 232 Th source. 
The energy spectrum of the events collected by CUORI- 
CINO can be seen in Figure [T] A more detailed descrip- 
tion of the experiment can be found in [Io| . 

CUORICINO 's geometry provides a unique opportu- 
nity to search for Oz/,3/3 and 2vf3[3 decay to the first + 
excited state in 130 Xe in an essentially background-free 
environment. This is due to the fact that these pro- 
cesses can be studied using a coincidence-based analysis 
by searching for two 7 lines of well defined energy. As 
can be seen in Figure [21 the decay to the first + excited 
state in 130 Xe differs from the one to the ground state 
in that it produces a gamma cascade. Given the Q-value 
of the decay, Q^/3=2527.5keV [TT| - fT3l ] . the two electrons 
are left with a total energy of 734.0 keV. The most prob- 
able de-excitation pattern, with a 86% branching ratio, 
proceeds through the emission of a 1257.41 keV and a 
536.09 keV gamma. Though 2z//3/3 and Oz^/3/3 decay both 
result in the emission of two electrons, the spectra of 
the sum energy of the two electrons differ drastically. In 
the first case, the two resulting betas have a continuous 
spectrum in the range (0-^734.0) keV, while in the sec- 
ond case, the result is just a monochromatic beta peak 
centered at 734.0 keV. Theoretical evaluations and ex- 
perimental limits for these two processes can be found 
in Table fl] It is important to note that the theoretical 
calculation for the half life of 2z//3/3-decay to the first ex- 
cited state + reported in Table Q] is not the one originally 



Decay Transition Theoretical (y) Experimental (y) 

~~ ~ 0+ ->• Of 7.5x10^ [17, ^8] > 3.1x10^ [lV] 

V _ + ^0 + (1.6-M5) xlO 23 [20] > 2.8xl0 24 [21] 

~ 0+ -5-0^ (5.1-M4)xl0 22 [15, 16£ > 2.3xl0 21 [16] 

V + ->0+ (1.7^70) xlO 19 [20] 7.0xl0 20 [22] 

a Corrected values for 15] (discussion in text) 

TABLE I. Theoretical evaluations (for mpp=\ eV) and exper- 
imental best limits (90% CL) for the half-life of 130 Te Or/fifi 
and 2i>pp decay. 



indicated in reference [15J since it was based on a wrong 
evaluation of the phase-space. The reported value is the 
one re-elaborated by A.S. Barabash [16| on the basis of 
the correct phase space factor. 



II. SEARCH STRATEGY AND EVENT 
SELECTION 

In this analysis, we consider only configurations in 
which the electrons are contained in the crystal where 
the decay takes place, and each de-excitation photon is 
completely absorbed in one crystal. With these require- 
ments, three different scenarios are possible (see Fig- 
ure [3]) . Scenario 1 takes place when both gammas es- 
cape from the original crystal. In scenario 2, the low- 
energy gamma (536.09 keV) is trapped in the original 
crystal with the betas, while the high-energy gamma 
(1257.41 keV) escapes. Scenario 3 is the opposite of sce- 
nario 2: the high-energy gamma (1257.41 keV) is trapped 
in the original crystal with the betas, while the low- 
energy one (536.09 keV) escapes. The signatures and the 
corresponding efficiencies are reported in Table [TTJ A fur- 
ther explanation of the calculation of the efficiencies can 
be found in Section UTT1 

The first-level analysis of the CUORICINO data is 
common to all physics processes to be studied and is de- 
scribed in detail in [21| . It starts from raw events and 
ends with a set of energy-calibrated hits associated with 
a time, a crystal, and other ancillary information, such as 
pulse shape parameters. In this phase of the analysis, a 
channel- and time-dependent energy threshold is applied 
to the data, based on the performance of each bolometer. 

For the processes studied in this paper, the analysis 
consists of defining signatures according to the three sce- 
narios reported in Table HH using them to select events 
from the CUORICINO data and evaluating the corre- 
sponding efficiencies from GEANT4-based Monte Carlo 
simulations [23[. 

Event selection criteria can be grouped into three 
categories: global, event-based and coincidence-based. 
Global and event-based cuts are not specific to this anal- 
ysis, and here we only outline them briefly (refer to [2"i| 
for details). Defined a priori, global cuts are used to dis- 
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energy (keV) 



FIG. 1. Single-hit (black line) and double- hit (red dashed line) energy spectra collected by CUORICINO in the range 
(500-^2700) keV. 



Decay mode Scenario 


Signature (energies in keV) 


MC 


Efficiency 
Instrumental 


Total 


0;/ 


1 

2 
3 


734 (/3) + 536 (7) + 1257 (7) 
1257 (7) + 1270 09 + 7) 
536 (7) + 1991 (p + 7) 


(0.60±0.02)% 
(2.29±0.04)% 
(1.41±0.03)% 


(86±2)% 
(90±1)% 
(90±1)% 


(0.44 ± 0.02)% 
(1.77 ± 0.04)% 
(1.09 ± 0.03)% 


2;/ 


1 

2 
3 


(0 -7- 734) 09) + 536 (7) + 1257 (7) (0.53±0.02)% 
(536 -r 1270) (/3 + 7) + 1257 (7) (3.04±0.04)% 
(1257 -T- 1991) (/3 + 7) + 536 (7) (1.28±0.03)% 


(86±2)% 
(90±1)% 
(90±1)% 


(0.39 ± 0.02)% 
(2.35 ± 0.04)% 
(0.99 ± 0.03)% 



TABLE II. Signatures and efficiencies for the three scenarios for Qv and 2v decay. We denote with the + sign the coincidence of 
energies released in different crystals. Efficiencies labeled as MC were computed based on Monte Carlo simulations. Instrumental 
efficiencies were computed based on CUORICINO data. Total efficiencies are given by the product of MC and instrumental 
efficiencies, times a factor of 0.86 to account for the branching ratio of the considered decay scheme (see Figure [3]). 

sion of non-physical pulses (electronic spikes or cryogenic- 
induced pulses) and physical pulses for which the energy 
is not estimated correctly (pile-up or excessive noise su- 
perimposed on the pulse). 

Coincidence-based cuts rely on the properties of a 
group of events that occurred within a fixed time window. 
Events can be selected based on the number of involved 
crystals, the spatial distance among them, the sum en- 
ergy or the energy of the single hits. In this paper, a 
100 ms time window was used to define coincident events. 
Physical coincidences induced by 130 Xe de-excitation oc- 
cur on much shorter time scales, but such a large time 
window must be chosen to account for the slow response 
of the bolometers. 

The coincidence-based event selection criteria were de- 
cided based on the scenarios described at the beginning of 
this section. Because the two electrons emitted in the 2v 
decay have a continuous spectrum, a wide energy window 
must be chosen for one of the crystals. This has the effect 
of introducing a much bigger background than is present 
in the analysis of the 0v decay mode. As a consequence, 
besides the criteria reported in Table lU additional re- 
strictions were applied to the events to be included in the 



6,38=2527.5 keV 



0; 



2? 



2; 



iS2_ 536.09 keV 



FIG. 2. Decay scheme for 130 Te, showing the energy levels 
(keV) and the branching ratios for the 7-rays [l4[ • 



card time windows in which one or more detectors per- 
formed poorly. This could happen because of external 
noise or cryogenic instabilities, which in turn result in a 
bad energy resolution. Event-based cuts allow the exclu- 



4 



(a) Scenario 1 



V 



(b) Scenario 2 




(c) Scenario 3 



FIG. 3. Possible capture scenarios. The blue lines represent 
the 1257.41 keV 7, while the red lines represent the 536.09 keV 
one. For each scenario, the available energy for the emitted 
/3s is 734.0 keV. 



2v analysis. To reduce random coincidences, a cut was 
imposed on the distance between the crystals involved 
in the events, as it was seen from the simulation that 
there is a low chance for the investigated processes to 
involve crystals that are far apart from each other. The 
most relevant background from physical processes is due 
to gamma rays that undergo a Compton interaction in 
one crystal and are then absorbed in another crystal. 
While the sum energy of these events is fixed, the energy 
released in each crystal has a continuous distribution. 
To reduce this background, events whose sum energy fell 
into a window of ±8 keV around the most intense gamma 
lines (1729.60 keV, 1764.49 keV, 1847.42 keV, 2118.5 keV, 
2204.21 keV and 2447.86 keV from 214 Bi, 2505 keV from 
60 Co and 2615 keV from 208 T1) were removed. 



III. ANALYSIS 

As stated in Section [III Monte Carlo simulations were 
used to calculate the efficiencies for the processes studied 
in this paper. This was achieved by comparing the num- 
ber of events passing the coincidence cuts to the total 
number of simulated events. The relatively low efficien- 
cies reported in Table [TT| arise from the fact that most of 
the gammas escape the crystals undetected and are ab- 
sorbed by inert materials surrounding them. Moreover, 
the signatures sought only consider the case of photons 
that are completely absorbed in one crystal, thus reject- 
ing events in which at least one photon is absorbed in 



one crystal after undergoing a Compton interaction in 
a different one. The computed values reported in the 
last column of Table [TT| also include inefficiencies due to 
event-based cuts, channel- and time-dependent energy 
thresholds, and discarded time windows in which one 
or more detectors were not performing properly (global 
cuts). Inefficiencies induced by channel based cuts were 
evaluated on the CUORICINO data in the same way dis- 
cussed in [2l| . The effect of global cuts was taken into ac- 
count by removing the simulated events lying in the time 
windows that were discarded from the real CUORICINO 
data, after said time windows were rescaled by the ra- 
tio between the total duration of the simulation and the 
real CUORICINO live time. The same procedure was 
used to associate energy thresholds to the simulation. 
Because the effect of global cuts was taken into consid- 
eration when determining the efficiencies, the exposure 
used in this work corresponds to the complete CUORI- 
CINO statistics without any subtractions: N( 130 Te)-t = 
9.5 x 10 25 y. 

Figure U shows the energy spectra obtained from the 
CUORICINO data after applying the event selection cuts 
described in Section [TTJ For each scenario, the spectrum 
was built as follows. Coincidence cuts were applied based 
on Table [Til requiring that the accepted events be in co- 
incidence with events satisfying each component hit of 
the signature except for the hit corresponding to the 
highest-energy 7. The signal search could then consist 
of a search in the resultant spectrum for evidence of the 
highest-energy 7 of the signature, which is the component 
with the lowest background. Moreover, the acceptance 
width for each cut was enlarged by ±10 keV with respect 
to the energies and energy ranges listed in Table [ill to 
account for the finite energy resolution of the detectors 
(a ~ 2 keV; see discussion below) . The energy windows 
used for the spectra were chosen to be much larger than 
the detector resolution, but small enough that at most 
one radioactive background peak was included, and the 
continuum could be assumed to be flat or linear. 



No evidence for a signal was found in any of the energy 
spectra. For the zero-neutrino decay mode, the back- 
ground is negligible, and no fit was performed. In this 
case, a condition of zero signal and zero background was 
assumed. In contrast, the background is not negligible for 
the 2v decay mode, and therefore a Bayesian maximum 
likelihood fit was performed for the 2v analyses. The 
best-fit curves are represented by the blue lines in Fig- 
ure [3J Depending on the scenario, different background 
models were adopted for the 2v spectra. The continuum 
was fitted with a constant (scenarios 1 and 2) or linear 
shape (scenario 3), while the possible additional peaks 
(1238 keV from 214 Bi for scenario 2, 511 keV for scenario 
3) were fitted with a Gaussian shape. The free parame- 
ters in the fit were as follows: the number of signal counts, 
the number of events from the flat background and the 
number of counts under the additional background peaks 
(scenarios 2 and 3). The energy resolution was fixed to 
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FIG. 4. CUORICINO energy spectra after the event selection cuts applied for the Of (left) and 2v (right) analyses. For the 
Of decay, the signal was expected at 1257.41 keV (plot A), 1270 keV (plot B) and 1991 keV (plot C) for scenarios 1, 2 and 3 
respectively. For the 2v decay, the signal was expected at 1257.41 keV for scenarios 1 (plot D) and 2 (plot E) and at 536.09 keV 
for scenario 3 (plot F). 



cr=1.8 keV. It was evaluated on the 511 keV peak and on 
the two 60 Co peaks at 1173 keV and 1332 keV that are 
visible in the CUORICINO energy spectrum (see Fig- 
ure [1} , and it was found to be comparable for all three 
peaks. A summary of the best-fit values for the 2i//3/3 
searches is reported in Table IIIII Systematic uncertain- 
ties were evaluated by repeating the fitting procedure 
with different background models, fitting ranges and en- 
ergy resolutions, and, compared to statistical uncertain- 



ties, they were found to be negligible. 



IV. RESULTS 

For each decay mode and for each of the three sce- 
narios, the posterior probability density function (p.d.f.) 
for the number of signal counts, P(JVg), was extracted 
using a Bayesian approach and assuming flat priors in 
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Scenario Ns 


N B 


[counts] 


[counts/keV] 


1 1.1±1.4±0.29 


0.12±0.03 


2 -0.4±6.6±2.6 


6.31±0.41 


3 -3.0±6.8±2.8 


6.73±0.33 



TABLE III. 2v analysis best-fit values for the number of signal 
(Ns) and background (Nb) counts. For Ns, both statistical 
and systematic uncertainties are reported. 



2 orders of magnitude, for both the Ov and 2v processes, 
with respect to the results of past experiments. It is 
worth noting that the new lower limit on the half life of 
the 2v decay mode is close to the upper bound of the 
theoretical calculation presented in Table |U A more clear 
picture will be available once CUORE, CUORICINO's 
successor [24|, comes online. This is due to the increase 
in target mass and improved background reduction that 
will be achieved in CUORE. 



the physical region (N s > 0). For the 0z/ decay mode, 
because there was no evidence of a signal and the back- 
ground was negligible, a Poisson p.d.f. for zero observed 
events was assumed for all three scenarios. For the 2v 
decay mode, the p.d.f.s were obtained as a result of the 
maximum likelihood fits on the spectra shown in Fig- 
ure |U For each decay mode, a global p.d.f. for the 
decay rate was obtained as the product of the three in- 
dividual p.d.f.s, Ptot(F) = YliPi(T)- In this formula 
Pi(T) = P l (N s ) ■ e l ■ N( 130 Te) ■ t, where the index i runs 
over the three scenarios and £j is the corresponding de- 
tection efficiency from Table |TTJ Systematic uncertainties 
were included in the Pi(T) according to the procedure 
described in [2l| . This resulted in the following half life 
lower limits: 

T 1/2 (2p/3/3*) > 1.3- 10 23 y, 90%C.L. 
T 1/2 (Qv(3(3*) > 9.4- 10 23 y, 90%C.L. 

These new limits represent an improvement of almost 
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